We present the Fourier parameter fit method, a new 
method for spectroscopically identifying stellar radial and 
non-radial pulsation modes based on the high-resolution 
time-series spectroscopy of absorption-line profiles. In con- 
trast to previous methods this one permits a quantifica- 
tion of the statistical significance of the computed solu- 
tions. The application of genetic algorithms in seeking so- 
lutions makes it possible to search through a large pa- 
rameter space. The mode identification is carried out by 
minimizing ^ ^ using the observed amplitude and phase 
\^ ' across the line profile and their modeled counterparts. 
, Computations of the theoretical line profiles are based on 
' a stellar displacement field, which is described as super- 
, position of spherical harmonics and that includes the first 
P5 • order effects of the Coriolis force. We made numerical tests 
^ ' of the method on a grid of different mono- and multi-mode 
• models for < £ < Ain order to explore its capabilities and 
limitations. Our results show that whereas the azimuthal 
order m can be unambiguously identified for low-order 
modes, the error of £ is in the range of ±1. The value of m 
^ ■ can be determined with higher precision than with other 
spectroscopic mode identification methods. Improved val- 
ues for the inclination can be obtained from the analysis 
of non-axisymmetric pulsation modes. The new method 
is ideally suited to intermediatley rotating S Scuti and j3 
Cephei stars. 
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1. Introduction 

Asteroseismology relies on a detailed comparison between 
the set of frequencies observed in a pulsating star and 
those predicted by theory. To do this successfully, we need 
to know the quantum numbers associated with each ob- 
served pulsation mode, and thus we need mode identifica- 
tions for as many frequencies as possible. The combined 
approach of spectroscopic identification methods, multi- 
color photometry, and frequency pattern recognition can 
identify both the harmonic degree £ and the azimuthal or- 
der m of a pulsation mode. The parameters of a stellar 
seismic model, such as the mass, the luminosity, the ro- 
tation rate, the effective temperature, the metallicity, and 
the age, can be tuned to obtain a fit between the theoret- 
ical and observed frequencies. In most cases the range of 
possible parameters is very large. By only considering the 
models that are consistent with mode identification, one 
is able to put a strong constraint on the internal structure 
of the star. For an extensive overview paper about photo- 
metric mode identification methods we refer the reader to 
Garrido (2000), and to Balona & Evers (1999) for a paper 
about the state-of-the-art approach of photometric mode 
identification. 

Profile variations of absorption lines formed in the pho- 
tosphere can be used to study the pulsational behavior of 
a star. The line-profile variations (LPV) are caused by the 
Doppler shift of absorption lines originating from differ- 
ent parts of the surface. High-resolution time-series spec- 
troscopy is a powerful tool for identifying the £ and m val- 
ues of a non-radially oscillating star by analyzing its LPV. 
In addition, these techniques allow determination of the 
inclination of the pulsation/rotation axis, the pulsational 
temperature variations, and the intrinsic oscillation am- 
plitude. If photometric amplitudes in different passbands 
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are also known, it is possible to calculate the ratio of flux 
perturbation to radial displacement, which is very sensi- 
tive to convection, thereby permitting selection between 
different convection models (Daszynska-Daszkiewicz et al. 
2003). 

In the past decade, considerable effort has been ex- 
pended to develop and test reliable spectroscopic mode 
identification methods. Basically, two kinds of approaches 
have been developed: one based on quantities integrated 
across the line profile and one based on the intensity vari- 
ations within each wavelength bin. The first approach is 
utilized in the moment method (Balona 1986, Aerts 1996, 
Briquet & Aerts 2003) and is applicable to the identifica- 
tion of low-degree modes {£ < 3). Different variants exist 
of the second approach, which is based on the Doppler 
image principle and requires sufficient rotational broad- 
ening of the profile. These include the intensity-period 
search (IPS) method (Schrijvers et al. 1997), where i and 
m are determined from the amplitude and phase distribu- 
tion across the line profile by applying empirically derived 
relations, and the pixel-by-pixel method (PPM) described 
by Mantegazza (2000). 

These spectroscopic mode identification methods have 
one thing in common: no statistical criterion that can 
quantify the significance of the obtained solutions can be 
calculated yet. For instance, the moment method (MM) 
quantifies the goodness-of-fit in elaborately weighted dis- 
criminants, which are not related to the observational un- 
certainties. Therefore it is often not possible to select be- 
tween different solutions if their discriminants have sim- 
ilar values. De Riddcr et al. (2003) report pessimistic 
prospects for defining such error estimates in the future 
because of the complex nature of the moments. 

We have developed a spectroscopic mode identification 
method, the Fourier parameter fit (FPF) method, which 
relies on a fit of the observational and modeled Fourier 
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parameters across the line profile. For every detected pul- 
sation frequency, the zero point, amplitude, and phase are 
computed for every wavelength bin across the profile by a 
multi-periodic least-square fit. The observational errors of 
these parameters can be calculated analytically from the 
least-square fit. This enables us to carry out a classical chi- 
square (x^)-test and to derive the statistical significance 
of our fits. 

The paper is constructed in the following way: wc first 
provide the formalism of the displacement field used for 
the computation of the LPV. The FPF method is pre- 
sented in the next section, followed by the description of 
numerical experiments applied for testing the method. 

2. Modeling of the line profile variations 

We assume that the displacement field of a pulsating star 

can be described by a superposition of spherical harmon- 
ics. Our description of the Lagrangian displacement field is 
valid in the limit of slow rotation taking the effects of the 
Coriolis force to the first order into account (Schrijvcrs 
et al. 1997). Since deviations from spherical symmetry 
due to centrifugal forces are ignored, our formalism is re- 
liable only for pulsation modes where the ratio of the ro- 
tation to the oscillation frequency 0/w < 0.5 (Aerts & 
Waelkens 1993). This limitation excludes realistic model- 
ing of rapidly rotating stars and low-frequency g-modes. 
For higher frequency p-modes, such as observed in many 
6 Scuti and (3 Cephei stars, the given criterion is fulfilled 
and a correct treatment is provided. 

The intrinsic line profile is assumed to be a Gaussian. 
This is a good approximation for shallow lines where 
the rotational broadening dominates other line broaden- 
ing mechanisms. A distorted profile is computed from a 
weighted summation of Dopplcr shifted profiles over the 
visible stellar surface. Additionally, wc take into account 
a parameterized variable equivalent width due to temper- 
ature and brightness variations across the stellar surface. 

Wc implemented the calculation of the pulsationally 
distorted absorption line profiles in a C-program code 
called LIPS, which can be obtained from the author upon 
request. The code is also capable to computing a displace- 
ment field whose symmetry axis is not aligned with the ro- 
tation axis, as observed in roAp stars (Kurtz 1982). Flux 
variations in different filters can be computed by consider- 
ing flux derivatives calculated from static atmospheres. A 
detailed description of the formalisms used for construct- 
ing the line profiles can be found in Appendix A. 

3. The xl of the FPF method 

Since the behavior of the LPV can be very complex due to 
the multi-periodicity of a star, the contribution of every 
single pulsation mode to the total profile variations must 
be separated. Thereby, it is possible to apply techniques 
developed for one- dimensional time-series data, if the in- 
tensity variations of every pixel or wavelength bin across 
the profile are analyzed separately. 



The search for pulsation frequencies can be carried out 
by applying a Fourier analysis for every pixel of the profile. 
For all detected frequencies it is possible, by means of 
least-square fitting, to determine the zero point Zx, the 
amplitude Ax, and the phase 0a across the line. We can 
thereby quantify the LPV caused by a pulsation mode by 
its Fourier parameters, which are a function of wavelength 
or Doppler velocity. 

Wc directly fit the observed Fourier parameters with 
theoretical values by applying an x^-test. The main dif- 
ferences to previous, similar mode identification methods 
are the following: 

— We utilize all available information on the Fourier pa- 
rameters, zero point, amplitude, and phase across the 
line profile. 

— The calculation of enables us to quantify the signif- 
icance of our fits. 

— The optimization is carried out with genetic algo- 
rithms, which permit the detection of local minima in 
a large multi-parameter space in much shorter compu- 
tational time than a grid allows. 

We calculate the theoretical line profiles with the pro- 
file generation code LIPS described in the first part of this 
paper. The theoretical Fourier parameters are computed 
from a least-square fit of typically ten synthetic mono- 
mode line profiles evenly sampled over one pulsation cy- 
cle. As for the observational Fourier parameters, here we 
also fix the frequency during the least-square fit of the 
wavelength bins and only derive Zx, Ax, and 

The observational variance cr^ of the intensity of ev- 
ery pixel across the profile must be derived very carefully, 
since the value of is very sensitive to this variance. 
Generally, is not constant along the profile, but in- 
creases towards its center due to the lower signal present 
there. The variance can be derived directly from the error- 
matrix of the multi-period least-square fit. 

Our approach is the following: we first determine the 
pulsationally independent stellar parameters usini, equiv- 
alent width Wo, and intrinsic line width a by fitting the 
observed zero point profile Zx- We assume Zx to be inde- 
pendent of the pulsation, which is true for a low ratio of 
the radial velocity amplitude to vsini. During the opti- 
mization process, we set the values of these parameters as 
variable in a narrow range set by their derived uncertain- 
ties. 

The parameters of the pulsation modes are derived 
only from the amplitude and phase across the line. The 
reduced xt is calculated from complex amplitudes in or- 
der to combine amplitude and phase information in the 
following way 
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where nx is the number of pixels across the profile, m the 
number of free parameters, A° and A* denote observa- 
tionally and theoretically determined values, respectively. 
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An = A\ COS (px and Aj — A\ sin 4>\ are the real and imag- 
inary parts of the complex amplitude, and a is the obser- 
vational error. Equationncan be easily modified if the ob- 
served amplitude and phase from photometric passbands 
is included for the calculation of xl- 

Since the amplitude and phase of a given wavelength 
bin are treated as independent variables, the variances are 
calculated from 

4,A = ^^(^a)' cos^ 0;, + a(0A)'^A sin' 0a (2) 

aj^^ = a{Axf sin' 0^ + ct(</)a)'A| cos' 0a. (3) 

In the case of a multi-periodic star, it is possible to 
optimize the sum of the xt values of all present modes by 
simultaneously taking common values for the pulsationally 
independent parameters inclination, vs'mi, a, and W. By 
doing so, the stellar inclination angle, and henceforth also 
the modes' intrinsic pulsation amplitudes, can be repro- 
duced better, since we take advantage of the fact that all 
observed modes are seen at the same aspect angle. Such 
an improvement in the identification of modes in multi- 
periodic stars has also been obtained by Briquet & Aerts 
(2003) for the case of the MM. 

We carry out the search for the best-fit model by means 
of a genetic optimization algorithm (Nissen 1997) by al- 
lowing the exploration of a large parameter space in much 
shorter computing time than can be done by grid calcu- 
lations. The basis of a genetic code is a stochastically se- 
lected set of models, which evolves in time towards better 
models by favoring those having a higher fitness (lower 
xt)- The strength of genetic algorithms is especially their 
ability to explore many local minima simultaneously and 
consequently to locate the absolute minimum of a complex 
multi-parameter optimization problem. 

4. Testing the FPF method 

We tested the functionality of the FPF method in a clas- 
sical way via numerical experiments by fitting simulated 
observations of mono-mode pulsations. We especially ex- 
amined to what precision different parameters of the input 
modes can be determined, with an emphasis on £, m, the 
relative displacement amplitude Os, and the inclination 
i (see Appendix A for a definition of these parameters). 
To check the stability of our solutions we made use of dia- 
grams, where x' is plotted as a function of every single free 
parameter. In such diagrams we mark the x'-values of all 
models computed during the optimization. Models having 
x' below a certain significance limit can be regarded as 
significant solutions representing the observations. 

In a series of four papers, Schrijvers et al. (1997), 
Telting & Schrijvers (1997a, 1997b), and Schrijvers & 
Telting (1999) have extensively tested the distribution of 
amplitude and phase across the line profile for different 
settings of mode input parameters. Therefore, we only fo- 
cus on the results of fitting the Fourier parameters and in 
this paper we only present the tests for low-degree modes. 



where temperature variations are not included in the mod- 
eling of the LPV and where the pulsational symmetry axis 
is aligned with the rotation axis. 

4.1. Low-degree modes 

We tested the quality of our fits for low-degree modes 
having < ^ < 4, all possible positive m- values with 
< m < £ and Ug = 0.001. Our synthetic line pro- 
files were computed by integration over a regular latitude- 
longitude surface grid consisting of 16200 visible elements, 
each with a size of 2°x 2°. For each test case we cal- 
culated irregularly spaced time series of 300 spectra for 
a timebase of 12 days with a pulsation frequency of 
/ = 12 d~^. The projected equatorial rotation velocity 
was set to vsini = 30 kms~^; furthermore, we set the 
intrinsic line width a ^ 12 kms~^ and the equivalent 
width W = 8 kms~^. These two values are based on ob- 
served values of the Fe i A5367.467 A line of the 6 Set 
star FG Vir. Other physical properties of the model (mass 
M — 1.85 Mq, radius R = 2.27 i?©, effective temperature 
Toft = 7250 K) also resemble those of FG Vir (Breger et 
al. 1999). We adopted the quadratic limb darkening coef- 
ficients Ua = 0.526 and ui, — 0.163 listed by Barban et al. 
(2003) linearly interpolated to the selected wavelength. 

To create conditions that are as realistic as possible, 
the raw synthetic spectra were convolved with a Gaussian 
instrumental profile resembling a spectrograph with a res- 
olution of R=60 000. Random noise was added to result 
in S/N=200. 

We derived the uncertainties a{Z\), a{A\), and (j{4>\) 
from the standard deviation of the residuals across the 
line profile after prewhitening the input frequency / and 
its first harmonic 2/. An estimate of usini=30.1 kms~^, 
a=11.9 kms"-'^, and W—8.0 kms~^was derived by fitting 
the zero point profile. Since the equivalent width W is very 
well constrained, we set it as a constant throughout the 
optimization, whereas the other two parameters were set 
as free in a range of ±0.5 kms^^ of the derived values. In 
Table n] the three best solutions from a fit to the Fourier 
parameters are listed for each tested mode setting. 

For selected modes we present x^-diagrams in Fig. ^ 
where the behavior of xt with respect to the free parame- 
ters is displayed. Every dot represents a model computed 
during the optimization. Some general conclusions, which 
apply to all tested pulsation geometries are the following: 



— For all tested pulsation geometries, the correct identifi- 
cation always had the lowest x' ■ There were sometimes 
ambiguities where more than one solution was below 
the 95 % confidence limit. 

— The value of m is in any case better constrained than 
£. The correct identification of m is always in first 
place, whereas some ambiguity often exists for typing 
£. The ambiguities in determinating £ and m increase 
for higher degree input modes. 



Zima : A new method for the spectroscopic identification of stellar non-radial pulsation modes 



5 



Table 1. Mode parameters derived from the application of the FPF method. The following parameter space was 
searched using genetic optimization (the values in brackets denote the range [min,max;A]): I [0,7;1], m [-7,7;1], 
a, [0.0001,0.005;0.0001],i [5°,90°;5°], ?;sini [29.5, 30.5; 0.1] kms"!, a [11.5, 12.5; 0.1] kms"!. For each tested displace- 
ment field, the three best solutions are shown. We used the following fixed model input values: pulsational amplitude 
as = 0.001, stellar inclination angle i = 40°, projected rotational velocity i)sini=30 kms~^, width of the intrinsic 
profile C7 = 12 kms~^, and equivalent width W = 8 kms~^. The 95 %-confidence limit of is 1.35. 
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— In most cases, the derived values for the amplitude and 
the inclination are correctly reproduced. For both val- 
ues a larger uncertainty exists for axisymmetric modes. 
There is a small dependence of xt with respect to the 
inclination for the £ = l,m = - mode, and any in- 
clination angle can be fitted satisfactorily. In contrast, 
the non-axisymmetric modes show a clear minimum of 
xi at or close to the input inclination. The best preci- 



sion for determining the inclination can be achieved for 
non-axisymmetric modes with £ = 1 or 2. Identification 
ambiguities of higher-degree modes increase the uncer- 
tainty. 



We now summarize the conclusions from our tests for dif- 
ferent kinds of pulsation modes. 
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Fig. 1. Application of the FPF method to synthetic LPV of three selected displacement fields (300 spectra, S/N=200). 
The left panels show, from top to bottom, the zero point, amplitude, and phase across the line of the input data (black 
crosses with error bars) and the two best fits (lines). The four panels on the right hand side show with respect 
to £ (top left), m (top right), the pulsational amplitude (bottom left), and the inclination angle (bottom left). Every 
single dot represents a model computed during the genetic optimization. The vertical solid line is at the position of 
the input parameter values and the horizontal line denotes the 95 %-confidence limit. 
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— Radial displacement: Both spherical mode param- 
eters, I and m, could be derived unambiguously, and, 
since for the radial mode no dependence of the ob- 
served amplitude on the inclination exists, the intrin- 
sic pulsation amplitude could also be derived correctly. 
When applying lower S/N values to the spectra, it is 
no longer possible to select between £ = and £ = 1, 
whereas only m = is still below the critical limit. 

— Axisymmetric displacement: The amplitude-phase 
pattern of these modes is very similar to that of radial 
pulsation, which often makes a distinction difficult. 
Again, the value of m can be determined precisely, 
but ambiguities exist in the typing of I and especially 
in the amplitude and the inclination. The last two pa- 
rameters directly depend on each other and a wrong 
inclination can be fitted by tuning the amplitude. 
Tesseral and sectoral displacements have in com- 
mon that the determination of the inclination and in- 
trinsic amplitude can be done with higher precision 
than for axisymmetric modes. Sectoral modes are more 
easily seen at higher inclinations, leading to a better 
identification due to the higher observed amplitude. 
An interesting aspect is that a misidentification with 
^derived = ^input + 1 Still generally leads to a good esti- 
mation of i. 

De Pauw et al. (i993) have carried out extensive 
numerical tests to explore the accuracy of the MM. 
Qualitatively, our conclusions for the FPF method agree 
very well with their results, but we were able to quantify 
them statistically for the FPF method. 

4.2. Dependence on S/N 

For two displacement fields, t = m = and ^ = 3, m = 2, 
we tested the effect of different S/N- values on the uncer- 
tainty of the mode identification. The stellar and pulsa- 
tional parameters were identical to those given in Sect. 14. ll 
and S/N-levels of 50, 100, 200, 300, 500 and 1000 were 
tested for a data set consisting of 300 spectra. 

As expected, higher S/N ratios lead to less ambigui- 
ties in the mode identification. The situation for the tested 
£ = 3, m = 2 is reported in Fig. |21 For the lowest adopted 
S/N levels, 50 and 100, respectively, no unique identifi- 
cation can be acquired for both i and m. The situation 
improves significantly for S/N levels of 200 and higher, 
where only the fit of the correct input I and m values has 
an xt value below the 95 %-confidence limit. These tests 
show that spectroscopic measurements dedicated to mode 
identification assuming a number of 300 spectra should 
at least have an S/N of 200 to be able to provide a well- 
constrained mode identification. 

The variances (t{A\Y and <j{(t)\)'^ of Eqs. El and 13 
are indirectly proportional to the S/N of the input data 
and also indirectly proportional to the square root of the 
number of measurements. Thus, doubling the S/N has 
the same effect on the mode discrimination as a fourfold 
number of spectra, i.e., more emphasis should be set on 
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Fig. 2. Minimum values of for all tested i and m val- 
ues of an input ^ = 3,m = 2 mode. Every tested S/N is 
represented by a different line. The solid horizontal line is 
at the 95 %-confidence limit. The calculations were made 
for a fixed number of 300 input spectra. 

high quality measurements than on the number of spectra. 
More important than a large data set is the time base of 
the measurements that determines the frequency resolu- 
tion, i.e., the ability to separate close frequencies detected 
in many 5 Scuti stars (Breger & Pamyatnykh 2006). 

4.3. Multi-mode pulsation 

The experimental results from the mono-mode profile op- 
timization with the Fourier parameter method are also 
applicable to multi-mode pulsators, if linear pulsation the- 
ory applies. The Fourier parameters can be determined by 
least-square fitting for every single frequency that can be 
detected from a Fourier analysis of the first moment or of 
the pixel-by-pixel variations in the line profile. Schrijvers 
et al. (1997) point out that additional harmonics, beat, 
and sum frequencies exist in the case of multi-mode pul- 
sation. In our numerical simulations described below, we 
find that their amplitudes are often too low for detection, 
especially in the case of 5 Scuti stars. This is also demon- 
strated for FG Vir by Zima et al. (2006), where among 15 
frequencies only one harmonic and one combination term 
were detected. These additional frequencies cannot be ne- 
glected for objects, such as RR Lyrae stars, showing higher 
amplitude and must be included in least-square fits to de- 
rive correct amplitude and phase values for the intrinsic 
frequencies. 

We tested our method on a five-frequency model con- 
sisting of low-degree modes. A similar test has already 
been successfully carried out for a comparison of the PBP 
method to the MM by Zima et al. (2004). 

To create conditions that are as realistic as possible, we 
generated 594 line profiles taking timings from the 2002 
FG Vir campaign (Zima et al. 2003), which had a time 
base of 80 days. We added noise (S/N « 200) consisting 
of the residuals after prewhitening of 15 detected frequen- 
cies and simulated the effect of phase smearing for an inte- 
gration time of 7 minutes. The resulting pulsation pattern 
of the five pulsation modes was calculated by summing 
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up the contribution of every single spherical harmonic. In 
these synthetic spectra we did not include temperature 
effects so as to speed up the process of optimization. 

We were able to detect all five input frequencies from 
a Fourier analysis of the first moment, as well as from 
the pixel-by-pixel variations, which are more sensitive to 
higher-degree modes. Due to the relatively high noise, no 
additional peaks related to combinations or harmonics of 
the detected frequencies could be found. The mode iden- 
tification of every single detected mode was carried out in 
the same way as described for the mono-mode profiles. 

Consistent with the mono-mode fits, again the az- 
imuthal order m can be identified with much higher preci- 
sion than £. For all modes xt of the correct identification 
of m is below the 95 %-confidence limit. The precision of 
the derived i is about ± 1 for the tested modes. 
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Fig. 3. Determination of the stellar inclination angle from 
the application of the FPF method to synthetic line pro- 
files (usini = 30 km s~^, i — 25°). The derived xt~ 
values are shown as a function of the inclination for differ- 
ent displacement fields (solid lines) and of a multi-mode 
fit to the non-axisymmetric modes (crosses), respectively. 
Axisymmetric modes are not well-suited to deriving the in- 
clination angle. A simultaneous fit to all non-axisymmetric 
modes provides the best estimate of the inclination. 

An important result is that the stellar inclination angle 
and henceforth also the modes' intrinsic pulsation ampli- 
tude is better constrained for the non-axisymmetric modes 
(see Fig.lSJ. For axisymmetric modes, the sensitivity of 
with respect to the inclination is very small due to the fact 
that we cannot decouple the intrinsic amplitude and the 
inclination and, since we do not know the intrinsic am- 
plitude, we can scale it such that any inclination can be 
fitted satisfactorily. Especially the fact that an axisym- 
metric mode might be a radial mode makes it necessary 
to judge inclination values derived from such a mode with 
great caution. 

An improved approach to deriving a more precise in- 
clination value can be carried out by a simultaneous fit to 



the Fourier parameters of all detected non-axisymmetric 
modes with common values of wsini, a, and W. By doing 
so, we take advantage of multi-periodic pulsation, assum- 
ing that all observed modes are seen at the same aspect 
angle. Figure |3 shows the comparison of the derived xt 
values for the five input modes as a function of the in- 
clination. Whereas no significant minimum exists for the 
two axisymmetric modes /i and /4, there is a clear de- 
pendence on xt and the inclination for the other modes. 
The simultaneous fit of the three non-axisymmetric modes 
yields an inclination value of 25 ± 10°, which is consistent 
with the input value. 

4.4. Comparison with the MM and the PPM 

We carried out a comparison of the FPF method with two 
other spectroscopic mode identification methods: the MM 
and the PPM. This comparison was accomplished by ana- 
lyzing the same synthetic data set with all three methods. 
We used the data of the comparison test described in Zima 
et al. (2004), which consists of a time series of synthetic 
line profiles distorted by a 6-frequency displacement field. 
The input modes have a low degree (0 < f < 3), and 
the profiles are computed for a star having a low v sin i of 
21 kms~^. Since the three methods quantify the quality 
of their fits, which we call deviation parameter, in mathe- 
matically different ways, the obtained values were normal- 
ized to the minimum derived value to permit a meaningful 
comparison. 

The striking difference between the methods lies in 
the sensitivity of the deviation parameters with respect 
to different solutions. This is illustrated in Fig. ^ where 
the computed deviation parameters of the three meth- 
ods are shown for one of the identified displacement fields 
{£ = m = 1). The minimum deviation value of the FPF 
method is more pronounced, especially for the azimuthal 
order m, than the minimum of the two other methods, 
thus implying an improved constraint on the mode iden- 
tification. The situation is similar for the other five iden- 
tified displacement fields. The value of m is better con- 
strained for every test case than the harmonic degree £. 

The higher sensitivity of the FPF method compared 
to the two other methods is mainly due to the fact that 
only the deviations from the mean profile are taken into 
account for the calculation of x^- Therefore, the fitting 
puts emphasis on the time variable parameters and conse- 
quently on the determination of £, m, a^, and i, whereas 
less emphasis is put on wsini, a, and W. 

We conclude that by means of the FPF method, es- 
pecially m can be determined with much better precision 
than with the MM and the PPM. The high sensitivity 
of with respect to the azimuthal order increases the 
probability that only one m-value lies below the 95 %- 
confidence limit. In our comparison we made the assump- 
tion that systematic errors in the continuum normaliza- 
tion do not occur, which is often not true for real data. 
In fact the FPF method and the PPM are much more af- 
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Fig. 4. Comparison of the identification of I and m derived from applying the FPF, PPM, and MM to synthetic line 
profiles. The best derived deviation parameters for each tested value of £ and m is shown for the different methods. 
For comparison, the deviation values have been normalized. The model input values t ~ m — 1 are marked by the 
arrows. 



fected by normalization uncertainties than the MM. Any 
such systematic uncertainties would reduce the sensitivity 
of the first two methods compared to the MM. The behav- 
ior of the deviation parameters with respect to i is quite 
similar for all three tested methods, leading to ambiguities 
in the identification of this parameter. 

5. Application limits 

We examined the limits for applying the FPF method in 
terms of projected rotational velocity and radial velocity 
amplitude. A major limit is set by the pulsation amplitude 
of the modes especially in a multi-periodically pulsating 
star. If the ratio of the projected radial velocity ampli- 
tude (RVA) and the equivalent width of at least one mode 
is too high, other modes are no longer projected linearly 
onto the Doppler velocity scale relative to the center of 
the stellar disk. Furthermore, a significant fraction of the 
power of the modes is transferred to their combination 
terms decreasing the power of the intrinsic frequencies. 
The Fourier values of such modes are distorted and can- 
not be represented by a single-mode model. The zero point 
profile is broadened and will lead to an overestimation of 
wsini. 

This limit applies especially to narrow-lined stars with 
relatively large amplitudes, such as RR Lyr (Kolenberg 
2002) or V Eri (De Ridder et al. 2004), where spectroscopic 
mode identification using Fourier values contradicts the 
results obtained from photometry. For these cases the MM 
is much better suited since it uses integrated flux across 
the line. 

We tested this problem for a four-frequency multi- 
mode model and calculated Fourier parameter diagrams 
for different intrinsic pulsation amplitudes of /i, keeping 
the amplitudes of the other modes constant. The resulting 
diagrams are reported in Fig. A significant reduction of 
the amplitude across the profile occurs for /2, /a, and 
for (RVA/usini)/j above 0.13. 

The average and zero point profile respectively are 
broadened by a large amplitude mode. The effect of the 
overestimation of vsiui due to the pulsationally caused 



Table 2. Overestimation of wsini due to pulsational 
broadening of the zero point profile by a radial displace- 
ment field. The subscript i denotes input values, whereas 
d denotes derived values. 



RVA/u sini 


[v sin i\ i 


0.003 


1.000 


0.016 


1.000 


0.032 


1.001 


0.065 


1.007 


0.100 


1.010 


0.163 


1.026 


0.327 


1.131 


0.653 


1.455 



broadening was explored by fitting the zero point profile 
with synthetic counterparts. The results are reported in 
Table 121 For different arbitrary values of RVA/u sini, the 
ratio of the input to the derived v sin i is given. Up to 
an RVA that is 10% of wsini, the error of the derived 
wsini is below 1%. For ratios above 30%, where an error 
larger than 10% is induced, it is better to derive v sin i by 
fitting the narrowest available line profile. 

Other mode identification methods that rely on the 
calculation of Fourier parameters, such as the PPM or 
the IPS method, are affected in the same way. For stars 
with large RVA relative to their rotational line broadening, 
these methods have to be applied with great caution. The 
MM is much better suited to the mode identification of 
such objects. 

6. Conclusions 

We have presented a new spectroscopic mode identifica- 
tion method, the FPF method, based on fitting the Fourier 
parameters of the intensity variations across the line pro- 
file. Our description of the displacement field is valid for 
a slowly rotating non-radially pulsating star and includes 
toroidal motion, i.e., first-order effects due to the Coriolis 
force. The implementation allows a modeling of pulsa- 
tional geometries where the pulsational symmetry axis is 
not aligned with the stellar rotation axis. Temperature 
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Fig. 5. Influence of a displacement field having a large radial velocity amplitude (RVA) on the amplitude and phase of 
other simultaneously excited pulsation modes. The data were calculated from an experimental four-frequency multi- 
mode model. The computations were made for different amplitudes of /i, whereas the amplitudes of the other modes 
were not modified. The abscissa represents the Doppler velocity normalized to wsini. 



variations of the stellar atmosphere and their impact on 
the equivalent width variations of a line can be taken into 
account. 

The strength of the FPF method is its ability to esti- 
mate the significance of the derived mode parameters by 
means of an test. For a given error threshold it is pos- 
sible to exclude a number of tested £ and m values for a 
pulsation mode. Furthermore, the stability of the derived 
fits can be determined by analyzing the trend of the 
values across the tested parameter ranges. Compared to 
previous methods, this permits a statistical quantification 
of the obtained solutions. We have also shown that the 
FPF method can determine m with better precision than 
other methods. This is especially important for improving 
the treatment of rotational effects in theoretical seismic 
models (Pamyatnykh 2003). 

The presented method is optimized for slowly to in- 
termediately rotating pulsating stars that have a large 
projected rotational velocity compared to the pulsational 
radial velocity amplitude. Such objects can be found es- 
pecially among 6 Scuti and f3 Cephei stars (Aerts & De 
Cat 2003). In a subsequent paper we demonstrate the suc- 
cessful application of the FPF method to high-dispersion 



time-series spectra of the non-radial pulsating S Scuti star 
FG Vir (HD 106384). 

Appendix A: Computation of the distorted line 
profile and flux variations 

A.l. The displacement and velocity field 

We assume an unperturbed model star to be spherically 
symmetric, in hydrostatic equilibrium, and in the absence 
of a magnetic field and rotation. The position of a mass 
element of such a star can be written in spherical coor- 
dinates (r, 0, (j)) defined by the radius r, the colatitude 
9 [0, tt], i.e., the angular distance from the pole, and the az- 
imuth angle </> [0, 27r]. Any shift of a mass element from its 
equilibrium position is given by the Lagrangian displace- 
ment vector ^ = {^r T S,e , £,4>) ■ This displacement modifies 
the initial pressure po: the density £»0i and the gravita- 
tional potential $0 as a function of r, 9, </>, and the time 
t. The linear, adiabatic perturbations of these parameters 
are governed by the four equations of hydrodynamics, i.e., 
Poisson's equation, the equation of motion, the equation 
of continuity, and the condition for adiabacity. 
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This set of differential equations is solved by assuming 
that r, 0, and <p depend on Yf{0, 0) e'"*, where Yp{0, 0) 
(Abramowitz & Stegun 1964) denotes the spherical har- 
monic of degree £ and of order m, w is the angular pul- 
sation frequency, and t the time. The spherical harmonic 
can be written as 



Yr{0,^) = Ni 



'Pj'"lfcos( 



imtj> 



(A.l) 



Here, ' denotes the associated Legendre function of 
degree i and azimuthal order m given by 

(-1)™ o.m 



and 



Nl' 



(-1)- 



'(2£+l) {t-\m\)\ 
47r {l+\m\)\ 



(A.2) 



(A.3) 



is a normalization constant. The definition of N'P changes 
from author to author, which must be taken into account 
when comparing derived amplitudes. 

We model uniform stellar rotation, including first- 
order corrections due to the Coriolis force, which gives 
rise to toroidal motion. The resulting displacement field 
in the case of a slowly rotating non-radially pulsating star 
cannot be described by a single spherical harmonic any- 
more. It consists of one spheroidal and two toroidal terms, 
which only have a horizontal component, and is given for 
an angular frequency uj in the stellar frame of reference 
and a time t by 




O0 sm ocp 
Id 9 ' 



sm0d(j)' 80 



(A.4) 



(Martens & Smeyers 
Schrijvers et al. 1997). 



1982, Aerts & Waelkens 1993, 
Note that the term proportional 



to Y^-^ is not defined for radial and sectoral modes. Here, 
Qsj denotes the amplitude of the spheroidal component of 
the displacement field, whereas at^£+i and at^i-i are the 
corresponding amplitudes of the toroidal components. We 
neglect the first order correction of the amplitude due 
to rotation, whereby the amplitudes of the toroidal terms 
can be approximated by the following relations 

ni_jm\ + i^_ 



at.i-i = a,sP 



UJ £+1 2£+l 
n£+\m\ 2 



2^-1-1 



l + (^+l)fc 



(A.5) 



The factor v in in Eq. I A. 41 is introduced in order to scale 
the normalization y/iTrN^ = 1, such that as represents 
the fractional radius variation for radial pulsation. 



The ratio of the horizontal to vertical amplitude, which 
allows the distinction between p- and g-modes, can be ap- 
proximated in the limit of slow rotation by the following 
relation (Schrijvers et al. 1997) 



— = ka + 2m — 
as LU 



l + fcp 
£{£ + 1) 



(A.6) 



where ah and Os are the horizontal and vertical amplitude 
and Ci^n is the Ledoux constant. Likewise, fep ~ 



denotes the ratio of the horizontal to vertical amplitude 
in the limit of no rotation. Here, M and R are the stellar 
mass and radius, and G is the gravitational constant. 

The displaced surface S = {Sr,Sg,S^) is calculated 
from the sum of the undisplaced coordinates of the unit 
sphere and the displacement vector by 



(A.7) 



Here, the temporal dependence of the displacement field 
is already incorporated in the definition of ^ from Ea. lA.4l 
If the pulsation axis is tilted by an angle of j with respect 
to the stellar rotation axis, the resulting displaced surface 
is shifted by a transformation of the spherical coordinates 

by 

cos 0j = COS 6 cos j + sin cos (j) sin j 




sin 0j = V / 1 — cos^ 6j 



cos (jjj = (sin cos (j) cos j — cos sin j) / sin 9j 
sin = sin (j) sin 9 / sin 0j . 



(A. 



In our implementation we divide a sphere of unit ra- 
dius into a grid of regular segments each having equal 
dimensions in A0 and Acj). The displaced coordinates of 
the surface S are evaluated for every single point of the 
grid. 

The rotation of the surface is implemented by moving 
the observer around the surface, which is fixed at the ori- 
gin of the coordinate system. At time t and the inclination 
angle i, the spherical coordinates of the observer then are 
O = (1, z, fit). By orthographic projection, we convert the 
three-dimensional surface coordinates to a plane normal 
to the line of sight to two-dimensional coordinates by 



Sr 



COS Sg 



sin Sg sin(S'0 — ilt) 
a i — sin Sg cos i cos{S^ — fit) 



(A.9) 

These coordinates are used to calculate the area A{0, 0, z) 
of the projected surface segments and for displaying the 
distorted surface on a screen. At each time step of the cal- 
culations, the visible part of the surface S is determined 
by computing the angle z between the line of sight and 
the numerically derived surface normal of every single seg- 
ment. Visibility is provided if < z < Deviations from 
spherical symmetry due to the surface displacement field 
are taken into account for the calculation of z. 
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We compute the velocity field of the displaced surface 
from the numerical temporal derivative of the surface co- 
ordinates. The advantage over an analytical approach is 
its higher flexibility if the description of the displacement 
field is modified, e.g., if additional rotational effects are 
incorporated. The velocity component V in the direction 
of the observer finally is given by 

V — Vx sini cos{ilt) + Vy smism{ilt) + Vz cos i, (A. 10) 

where {Vx, Vy,Vz) are the Cartesian coordinates of the de- 
rived surface velocity field. 

A. 2. Construction of the distorted line profile 

We assume that the intrinsic line profile is a Gaussian, 
which may undergo equivalent width changes due to tem- 
perature variations. The distorted line profile is calculated 
from an integration of an intrinsic profile over the whole 
visible stellar surface, which - for computational purposes 
- numerically results in a weighted summation over the 
surface grid. We define the intrinsic Gaussian profile in 
a surface point having the line-of-sight velocity V (see 
Ea. rOT)|l as 



/(t;,Teff,logg) 



1 



6F 



1 



W^int(Teff) 



(A.ll) 

Here, v is the velocity across the line profile, 6F/F, which 
will be defined in Section IA.3I takes the surface fiux of 
the emitting segment into account, Wint{Tcs) is the equiv- 
alent width as a function of the effective temperature (see 
Ea. lA.13ll : and a is the width of the intrinsic profile. Note 
that in our definition the value of a is higher by a factor 
of compared to the definition by Aerts et al. (1992) 
generally adopted for the MM. 

The distorted profile 3 can be calculated by sum- 
mation over all visible segments on the surface grid of 
{6, (p) weighted over the projected surface and adopting a 
quadratic limb darkening law by 



E 



5{e, </., z)A{e, 0, z) 1 - K(i + ub{i - ^^y 



^/(ti. Toff, log 5) 



(A.12) 



Here, S{9, (p, z) — 1 if the surface element at {9, (p) is visible 
and 6{9, (p, z) = otherwise. Also, ^ = cos z is the cosine of 
the angle between the surface normal and the line of sight, 
Ua and Mb are the first and second order limb darkening 
coefficients, and A{d, cp, z) is the projected surface area of 
the surface element. 

The response of a line's equivalent width to local tem- 
perature changes is dependent on the involved element, its 
excitation, and the temperature in the zone where the line 
originates. In order to take this effect into account, a vari- 
able equivalent width of the intrinsic line profiles must be 



considered for calculating the distorted profile. Since there 
is no phase shift between 5Wint{T) and 5T we can write, 
following Schrijvers & Telling (1999), 



Wint(reff) = Wq[1 + awST.s), 



(A.13) 



where aw is a parameter denoting the equivalent width's 
linear dependence on STcs, which can be approximated 
for STcs ^ 1. A description of the temperature variations 
will be provided in the next section. 

A. 3. Flux variations 

For calculating the local temperature, surface gravity, and 
flux variations, we closely followed Balona (2000) and 
Daszynska-Daszkiewicz et al. (2002). Since the fiux vari- 
ation SF/F is mainly a function of Tcs and log 5, we can 
write in the limit of linear pulsation theory 



ST, 



OiT- 

6R 



cff 



an 



arf-e 



(A.14) 



where ar and ag given by 



ax 



dlogF 



dlogF 

dlogTcsJg " \d\oggyj.^^ 



and a a 



(A.15) 



are partial derivatives of the flux, which can be calculated 
from static model atmospheres for different passbands. 

Here, i?o is the unperturbed radius, G denotes the 
gravitational constant, < p > is the mean density of the 
star, / the absolute value of the complex fn + ifi, and ipf 
the phase lag of the displacement between the radius and 
temperature eigenfunctions. Then / describes the ratio of 
flux to radius variations, which can be transformed into 
the ratio of temperature to radius variations due to the 
fact that the flux is proportional to T"*. 

The total flux ^ is derived by weighted summation over 
the visible surface 



d{9, cp, z)A{9, 0, z) 1 - [ua{l - m) + "b(l - M)' 



SF_\ 



(A.16) 



where the symbols have the same meaning as in Eq. IA.12I 
By deflnition, / reflects the temperature change due 
to the radius variations during pulsation. Then dT^s in 
turn governs the equivalent width variations of a line and 
the flux variations. Therefore, we can determine empirical 
values of / by two independent approaches. The derived 
values of / yield information about sub-photospheric lay- 
ers where the thermal time scale is on the order of the 
pulsation period. These layers are only poorly probed by 
seismic modeling of only frequency values. 
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